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Abstract 

The free propagator for the scalar A0 4 -theory is calculated exactly up to 
the second derivative of a background field. Using this propagator I com- 
pute the one-loop effective action, which then contains all powers of the 
field but with at most two derivatives acting on each field. The standard 
derivative expansion, which only has a finite number of derivatives in each 
term, breaks down for small fields when the mass is zero, while the ex- 
pression obtained here has a well-defined expansion in 0. In this way the 
resummation of derivatives cures the naive IR divergence. The extension 
to finite temperature is also discussed. 



E-mail: per.elmfors@cern.ch 



1 Introduction 



There are a number of methods for computing approximations to the effective action 
depending on which parameters can be considered as small. If the field amplitude is 
small it may be enough to compute a small number of n-point functions, with the 
advantage that they can be computed to all orders in the derivatives. For instance, the 
two-point function can be calculated to all orders in momentum, yielding corrections 
to the dispersion relation. For configurations with large fields all n-point functions 
become important and have to be resummed. This we can do for a constant field, thus 
assuming all derivatives to be zero, and we get the effective potential. Even though 
the effective potential is calculated for a strictly constant field, the effective action, 
expanded to the zeroth order in derivatives, is obtained from it by simply replacing the 
constant by a slowly varying x-dependent field. By doing so it is assumed that the field 
varies so slowly that it can locally be approximated by a constant. Effects that have to 
do with global properties, such as phase separation in spontaneously broken theories or 
instanton contributions, are completely neglected in the quantum corrections, although 
there is no direct assumption that the amplitude of the field variation must be small. 
In a derivative expansion of the effective action the next step is to calculate the wave 
function renormalization to all orders in the field without derivatives. The reason why 
the effective potential and the wave function renormalization can be calculated to all 
orders in the field from a one-loop calculation is that the propagator can be constructed 
exactly in the presence of a constant background field. It is, therefore, interesting to 
see if one can go further in the local approximation while still being able to find the 
exact propagator. In this paper I use Schwinger's proper-time method Q to compute 
the propagator exactly up to the second derivative of the field and use this in turn to 
compute the effective action and the effective equations of motion. 

I restrict the calculations here to a toy model, the A0 4 -model, to study general 
properties of the method. It would be interesting to extend the method to more 
realistic theories, in particular gauge theories. The effective potential has been used 
to a large extent in the electroweak model to calculate bubble properties at the phase 
transition. The wave function renormalization constant Z((p) has also been computed 
[0]. There are two main problems with that procedure. First, it is very sensitive to 
the choice of gauge fixing parameter, because the usual effective action suffers from 
gauge dependence away from its stationary points. (The Vilkovisky-DeWitt effective 
action might be the correct way to get around this problem.) Secondly, Z(<p) diverges 
when the field goes to zero because the transverse gauge field is unscreened at high 
temperature. In this paper I discuss a way to cure the IR sensitivity by including 
higher derivative terms. The idea is that if the field is much smaller than the scale 
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set by the derivative, then the derivative expansion does no longer hold and higher 
derivatives have to be included exactly. On the other hand, when the derivatives are 
non-negligible, the field cannot be small in any sizeable volume and in this way the 
derivatives regulate the IR sensitivity. The approximation is still local and there is 
no hope that it would cure the imaginary parts that occur in a spontaneously broken 
theory when the mass squared is negative. It may, therefore, be difficult to apply the 
method directly to bubble wall calculations in the electroweak theory. 



2 Local expansion of the effective action 

The basic equation to use when calculating the effective action is the generalized tad- 
pole equation || 

™ = r«, (!) 

which gives the relation between the effective action T[<f)] and the one-point func- 
tion rW[0;x] computed in the background of <p(y). To find one has to inte- 
grate rW[0;x] with respect to <fi(x), but the equation of motion is given directly by 
T^ 1 ' [(/); x] = 0. The Feynman rules needed to compute the one-point function are 
derived from the Lagrangian in the non-trivial background 

m = \ (W) 2 - mV) - j y <f , (2) 



The tadpole equation to one loop is then written as 

ir u, [M = ^ + i HA ^. ))W __l_ w , (4) 

with M 2 (0) = m 2 + \<j) 2 (x)/2. This equation can be integrated to give 

T[(/>] = T d + l -Tr x ln(p 2 — M 2 + ie) . (5) 

Because of the trace over one-particle states labelled by x, there is no ordering problem 
when doing the integration, even though p and 4>{x) do not commute. From Eq. (|5|) 
we identify the one-loop Lagrangian and use the proper-time representation to write 
it as 

C^\x) = -i r -(x\e is ^- M2+i ^\x) = - l - [°° -(x(s)\x(0)} . (6) 
2 jo s 2 Jo s 
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This expression has, of course, the usual UV infinities which we have to renormalize at 
the end. 

In a general background it is not possible to calculate (x\ exp[is(p 2 — M 2 + ze)] |x) 
exactly since [p 2 , M 2 (0(x))] 7^ 0. If we think of Eq. (^) as a purely quantum mechanical 
problem, the non-trivial part of the amplitude from \x) to (x\ comes from quantum 
fluctuations. That is, the virtual particle that propagates from x back to x probes 
a certain neighbourghood of x related to the Compton wavelength. Expanding the 
background field, i.e. M 2 (0), in derivatives should lead to a good approximation if 
these fluctuations are not too large. Notice that this is different from a derivative 
expansion of the amplitude itself, which is a complicated function of M 2 . To second 
order in derivatives we define 



CO 



d^dyM 2 = Xd^ducj) + X^d^c 



(7) 



When computing the Lagrangian at a point x we expand the position operator in 
M 2 (0) around that point (ihj + i, x\x) = 0), and define the corresponding bilinear 
Hamiltonian (suppressing the +ie term in the sequel) 

M 2 {<p{x + x)) ~ M 2 + x^F^ + -Fw^f , 



H 



p 2 — M 2 — xF 



-XUJ X 



(8) 



using a matrix notation in the second line of Eq. (^). Following Schwinger [jl] the 
problem has now been reduced to solving the differential equation 



±(x"(s)\x'(0)) = i(x"(s)\H\x'(0)) , 
(x"(0)|x'(0)) = 5(x"-x') . 



(9) 



In order to write Eq. (g) as an explicit differential equation, H(p(0), x(0)) has to be 
rewritten as a function of x(s) and x(0) using the equations of motion for p and x 




H(p(s),x(s)), 



x(s) 
p(s) 




■ (10) 



The solution is 



x(s) +uj- l F 



cosh(sv2w) sinh(sv2u?) 
sinh(s-\/2u;) cosh(sy / 2ct;) 



x(0) +oj- 1 F 

'Sp(o) 
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and the differential equation 
d 



ds 



[x 



F" 



1 



2a; sinh 2 (s 

isM 2 ((f)(x')) 




is explicitly written as 



F" + F' 



1 cosh (syluj) 



2uj sinh 2 (s 
a7cosh(sv / 2cJ) \ 



F' 



F „ cosh(sV2aT) p , 



2 sinh(sv2a;) / 



2a;) a;sinh [s 

x"(s)\x'(0)) , 



2a;)" 



(12) 



where we have expanded M 2 (0) around x' . This choice of expansion point will not 



matter when we later take x' 



x 



x. The notation F" and F 1 means F„ 



x 



and 



F^(x'). The non-trivial functions of u should be defined by power series expansions 
in o; M v with suitable contractions with F^. The solution of Eq. ([12]) in <i-dimensional 
Minkowski space is 

1/2 

(x"(s)\x'(0)) = e-' l ^ d - 2) (47is)- d/2 det 



exp 



Finally, taking x' 
C w (x) 



x 



+ F" 
x we find 

dS /A \ 



^sinh(s\/2c^) / 
^coth^v 7 ^)^, F , 

1 



coth(sv2w) 



{2uofl 2 
1 



(2o/ 



3/2 



2c<j c<j sinh(s\/2a;) 



F' - sM 2 



(13) 



exp 



-is 



M 



-rf/2 e -if( 



2cj 



- 2 ) det 



1 - 



2a; 



1/2 



v sinh(s\/2a;) / 
tanh 



(14) 



Notice that the integrand is well-defined as a power series in o; M v even when uo^ v is 
not invertible. The integrand has possible poles on the real and imaginary axes, but 
only there since eigenvalues of u are real. We now want to deform the s-contour to 
the negative imaginary axis. This is possible if the arc at infinity gives a vanishing 
contribution. In our case we must require M 2 > |F M (a;~ 1 ) A1 ^F / . If we assume the 
eigenvalues of uj^ v to be positive, the condition on M 2 is the same as saying that the 
mass squared, in the quadratic approximation M 2 (x + £) ~ M 2 (x) — F M £ M + |£ M o; Ml , 
is positive everywhere, not only at x. We now assume this condition to be satisfied 
and perform the contour deformation. If the above condition is not satisfied, we would 
expect an imaginary part since the particles can then (in the quadratic approximation) 
propagate to a point where they become tachyonic. Such an imaginary part would be 
an artefact of the approximation rather than a physical reality if the full M 2 (0) indeed 
is positive everywhere. 
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In 4 dimensions (d = 4) Eq. (|i~4"D is infinite as it stands and has to be renormalized. 
This can be done in a standard manner by adding and subtracting 



32vr : 



;(V)" 



00 







e 3— e 



1 + 



-sAI 2 



(15) 



It is not necessary to renormalize the term proportional to u a a since it is finite after 
a partial integration (the divergent part is a total derivative). We add it anyway to 
make C^\x) finite before any partial integrations. The regularized expression in 4 — 2e 
dimensions is then 

1/2 

x) = / — det 



327T 2 



x exp < sF — 1 
P I 2u\ 





-sM 



Anfi 2 J 32tt 2 ~ a '"4vr/i 2 ' ^ 

Even though Eq. (|16|) looks manifestly real, there might be some imaginary contribu- 
tions from poles on the real s-axis, depending on the eigenvalues of u. The contour 
should go slightly above those poles. 

At this point it is interesting to start to series-expand Eq. (|T6| ) in derivatives to see 
if it coincides with known results. We then find before any partial integration: 




to 



^ln 



M 2 



1 F»F„ 



32vr 2 6 471-^ 32tt 2 M 2 



F^FV, 
72M 6 



+ 



+ 



144M 4 180M 4 



(17) 



The effective potential (no derivatives) and the wave function renormalization (two 
derivatives) come out correctly up to partial integration. However, comparing with 
e.g. or |J one finds an apparent discrepancy for the four-derivative term (in [|| 
there is also an overall sign error). The difference is a total derivative up to a term 
which contains a field with three derivatives: 




M 6 



80M 4 / 180M 4 



(18) 



In the present approximation we neglect the last term in Eq. (|18"D so that the result 
here is actually consistent with [|| [|. This observation reminds us that the local 
effective Lagrangian density in the effective action is not uniquely defined since we can 
add total derivatives without changing the equations of motion. Only the integral over 
space-time or the equations of motion have intrinsic meaning. 
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3 Zero field limit for m 2 = 



At the end of Section |2| we saw what a typical higher order term in the derivative 
expansion looks like (see Eq. (|17|)). It is easy to see that, for dimensional reasons, the 
higher the power of derivatives is, the higher the power of M 2 in the denominator is. For 
a massless theory, M 2 = A0 2 /2, so for small this expansion obviously breaks down. 
In the derivative expansion, corrections to derivative terms are computed assuming 
that the field is constant. Since the corrections diverges we have to resum the series. 
After all, if d^cj) is non-zero then cannot be zero everywhere. The propagator we 
calculated in Section [2] contains the resummation necessary to make the 0^0 limit 
meaningful. It does not really make sense to talk about = in the effective action 
since for instance (9^0(9^0 can be partially integrated to 0(9^*9^0, which would then 
be zero although it does contribute to the equation of motion. In order to avoid the 
problems of possible partial integration and mixing of derivatives we shall study the 
equations of motion directly. From Eq. (^) we have the unrenormalized equation of 
motion 

Ij- ~ 2 So = ° • (19) 

It may at first look strange that the whole one-loop contribution in Eq. (|T9| ) (using 
Eq. fll3|) ) comes from varying Eq. @ (using Eqs. (H0)) on ly with respect to the in 
M 2 and not in F or u. But, rewriting Eq. (|5|) as 

VM - - / -e-^W / d 4 y (y\ exp[ts(p 2 - M 2 (y) + M 2 {x))]\y) , (20) 
2 Jo s J 

one sees that the linear variation with respect to only gets contributions from the first 
factor e~ lsAl2 and not from the operator part in M 2 (y). Partial integrations in Eq. (|5|) 
would, of course, spoil this property. Also the truncation of including only the second 
derivative, as in Eq. flliD , spoils this property in general. We can verify explicitly from 
Eq. (|TT| ) that up to two derivatives 

5T dC „ dC „ „ dC dC 



5M 2 DM 2 ^dF^^du^ dM 2 ' (21) 
so that is does not matter whether it is only M 2 that is varied with respect to 



or also F and u. In the case of four or more derivatives Eq. ( 21[) no longer holds 
since for instance a term like d A M 2 /M 2 would contribute to the equations of motion 
in the present approximation, but it is not included in the effective Lagrangian in 



Eq. (14). This shows again that one should discuss the equations of motion where 



these ambiguities do not occur. 



6 



To simplify the analysis we shall now assume that the background field is static and 
that the eigenvalues of Uij are positive. This is natural since M 2 = X(f) 2 /2 is positive 
for real and has to be growing (or constant) in all directions away from a point where 
= 0. We continue to suppress the space-time indices but in this section they should 
be understood as only space indices with Euclidean metric. Thus, we define F and u 
by Fi and Uij, with all indices as subscripts, so there is a sign difference with lo^ v . In 
order to do a series expansion of Eq. ([0J) in powers of we write the one-loop part, 
after renormalization, as 



rW(x) = - 



sinh(sv 2uS) 




1/2 













exp 









tanh( 




(22) 



The two terms in the square brackets are separately finite for small s, and we shall 
treat them separately. The first one can actually be expanded directly in since the 
determinant makes the integral convergent for large s. The second term has to be 
treated with more care. One way to evaluate it is to divide the s-integral into two in- 
tervals, [0, A] and [A, oo], where A0 2 <C A -1 <C H^H 1 / 2 . Doing suitable approximations 
and partial integrations on each interval, we finally obtain 



rd), 



X 



A0 



32vr 2 




2 , W2ca 
-j= tanh(- 



(a 2 2 r-g{s)d+ 
I Jo s 

■ r*(i-g(s))+tf (--+ r^in(4v 

Jo s 2 \ e Jo 




dg{s ) 
ds 



where 



g(s) = det 



s\/2uj 



1/2 



,(23) 



(24) 



v sinh(sv2a;) / 

It is interesting to see that Eq. ( p3|) has a well-defined and finite expansion in powers 
of 0. The reason why it is IR-finite, i.e. why it is convergent for large s, is that 
g(s) — > rapidly. We do not see immediately from the treatment here how higher 
order derivatives of M 2 would affect the result, but we can expect the correction to be 
finite on physical grounds if there are no instabilities. 
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4 Finite temperature effective action 



There are two main formalisms for doing finite temperature field theory, the real-time 
and the imaginary-time formalisms, and both can be used to generalize the calculations 
in Section 0. We shall discuss both formalisms here and use them for different purposes. 



Real-time formalism 



The generalization to the real-time formalism is done by using the thermal propagator 
when calculating the tadpole in Eq. (^). Although a doubling of the number of degrees 
of freedom is generally needed for a consistent loop expansion, only the (ll)-part of the 
propagator is needed to calculate the tadpole at one loop. The thermal (ll)-propagator 
is given by 



+ 2f B (p )Re- 



(25) 



p 2 -M 2 + ie Jy ^ J p 2 -M 2 + ie' 

The real part can be represented by extending the s-integral in Eq. (||) from — oo to 
oo, but the phase of the normalization in Eq. flllf ) is only valid for s > 0. Instead of 
keeping track of the phase factor we continue with explicitly taking the real part. The 
order of /s(po) and M 2 (0(x)) in Eq. ( p5|) is not well-defined, unless the background 
field is time-independent. It also makes sense to limit the considerations to such 
fields since the system will not remain in equilibrium otherwise. It is conceivable to 
treat small time-dependent perturbations of an equilibrium background, assuming the 
perturbation to be small, but here we also want to deal with large field amplitudes. 
Thus, assuming F = and oj^q = 0, we get 



-Re 



2tt 



ds 



(4ns) 



-3/2 



X 



x exp 
1 



2u 



+ 



dot i 77_ ) 

\sin(sv2c<;) / 
isipl - M 2 ) 

dp fB(po)(p l 



1/2 



exp 




sy 2uj 



tan 




12tt 5 



M 2 ) 3 / 2 Q(p 2 - M 2 



(26) 



The last term in Eq. (|26|) is the effective potential, which comes out when we regularize 
for small s and use a T-function to define a finite value. 

This formal expression of the thermal effective action is rather difficult to use in 
general, since we would like to deform the s-contour to the positive imaginary axis for 
large po, but then we encounter all the poles on the positive real axis (the s-contour as 
it stands goes slightly below the real axis). At least if we assume the eigenvalues of 
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to be positive, which we want for stability. The situation is reminiscent of electrons in 
a constant magnetic field where the poles correspond to Landau levels, while here 
they are related to oscillations in a harmonic oscillator potential. 

As in Section ^ we want to compare Eq. fl2"E| ) with known derivative expansions. 
The effective potential is easily recognized to be correct from the last term in Eq. (26). 
The wave function renormalization comes from two terms, u>u and i^-Fj. A naive 
expansion in FjFj of the exponent in Eq. ( |26| ) leads to a divergent integral, but if one 
first makes a partial integration with respect to po it becomes finite. In the end the 
total contribution agrees with the high temperature expansion in ||. 



Imaginary-time formalism 

The transition to the imaginary-time formalism is most easily performed in a po rather 
than xo representation of the amplitude (x"(s)\x'(0)). In the limit of coinciding points 
we can write 

(x{s)\x(0)) = f ^dx' o (x[ ) ( S ),a; J ( S )|xo(0),x(0) J )e^^-o) j (27) 

J 271 

and replace / dpo by VItxT J2 n {Po i2nTn), as usual. Let us now study the equation 
of motion and see what happens in the <fi —>■ limit when m 2 = 0.0 The full one-loop 
correction to the equation of motion is 



r«[0] +r« [<f>\ = Hdse-^T V (4vr S )- 3 / 2 det 

2 Jo 



x exp 



- ls \^2nTn) 2 + M 2 -F^^l-^=t a n( S -^))F) . (28) 




We shall only study the leading contribution for small 0, so we can put M 2 = F = 0. 
The IR convergence is anyway governed by u. The zero temperature part of Eq. (|2~8|) 
can be extracted after a Poisson resummation, and in the remaining piece we deform 
the contour by taking s — > —is. Thus the finite value of VP for small is 

To get some more concrete information out of Eq. (p9| ) we can take the high temper- 
ature limit where the determinant can be approximated by 1. As expected we find 

2 To be consistent in the high temperature limit we would have to resum the thermal mass correction 
AT 2 /24. A zero effective mass is found at a second order phase transition in a spontaneously broken 
theory where the resummed mass is m 2 (T) = —fj, 2 + AT 2 /24, /i 2 > 0. This is really the situation we 
imagine when we write m 2 = here. 
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rl 1 ^ = —\T 2 <p/24, which simply is the thermal mass shift. Therefore, we also want the 
next to leading term in order to get something new. Again, a straightforward expansion 
of the determinant does not work since it gives divergent remaining integrals. This is 
not surprising since it is the large-s behaviour of the determinant that should make 
the expansion IR-convergent. One way to get around the problem is to first subtract 
the 1 from the determinant and then use the formula 



n=l 



-n 2 /s _ ±, 



-l) + v /^]Te-^, 



(30) 



n=l 



and a transformation s — > 1/s to restrict the integration interval to [l,oo]. We then 
obtain for the subleading contribution in T 



AT 2 



24 



\4>T 2 

4:71 Jl 



ds 



A0T 2 
' 8tt 2 



fds (g( 



-47tT 2 s' 

S s 



,-3/2 



•4ttT 2 



- 1 



n=l 



AttT 2 ' 



1 



-3/2 



using the notation in Eq. (U%- The first term on the right-hand side of Eq 
to zero in the high temperature limit, but the second term finally gives: 



r 



Mi 



AT 2 
~2T 



AT 
16^72"' 



ds 



s 



3/2 



det 



sy2u 



1/2 



K sinh(sv2a;) 



1 



(31) 
goes 

(32) 



We can then conclude that also the finite temperature part of the equation of motion 
is well-behaved when <ft — > as long as u has positive eigenvalues. 
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